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ABSTRACT 

We report on the discovery by the Swift Gamma-Ray Burst Explorer 
of the eighth known transient accretion-powered millisecond pulsar, SWIFT 
J1756.9— 2508, as part of routine observations with the Swift Burst Alert Tele- 
scope hard X-ray transient monitor. The pulsar was subsequently observed by 
both the X-Ray Telescope on Swift and the Rossi X-Ray Timing Explorer Pro- 
portional Counter Array. It has a spin frequency of 182 Hz (5.5 ms) and an 
orbital period of 54.7 minutes. The minimum companion mass is between 0.0067 
and 0.0086 M ? , depending on the mass of the neutron star, and the upper limit 
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on the mass is 0.030 M- (95% confidence level). Such a low mass is inconsistent 
with brown dwarf models, and comparison with white dwarf models suggests that 
the companion is a He-dominated donor whose thermal cooling has been at least 
modestly slowed by irradiation from the accretion flux. No X-ray bursts, dips, 
eclipses or quasi-periodic oscillations were detected. The current outburst lasted 
~ 13 days and no earlier outbursts were found in archival data. 

Subject headings: binaries: close — pulsars: general — stars: neutron — white 
dwarfs — X-rays:binaries 


1. Introduction 

It is now well established that a subset of neutron stars in low-mass X-ray binary systems 
are millisecond pulsars powered by accretion from either a brown dwarf or degenerate com- 
panion in a very close binary configuration. Such systems are believed to be progenitors of 
millisecond radio pulsars (Alpar et al. 1982). Before the discovery of SWIFT J1756.9— 2508, 
there were seven known accretion-powered millisecond pulsars (see Kaaret et al. 2006, for a 
full listing). The sample is now becoming large enough to start to study classes of these ob- 
jects which can in turn lead to a better understanding of the properties of neutron star X-ray 
binaries and their relationship to millisecond radio pulsars. The known accretion-powered 
millisecond pulsars all appear to be in very close orbits (with orbital periods ranging between 
41.1 and 256.5 min). For all but one such system, XTE J1814— 338 (Bhattacharyya et al. 
2005), orbital constraints rule out a hydrogen-burning main sequence companion, and for 
four systems, even a brown dwarf companion is not consistent with the orbital parameters. 
The application of recent models of white dwarf donors in ultracompact binaries (e.g. Deloye 
k Bildsten 2003; Deloye et al. 2007), shows that, within certain temperature ranges, white 
dwarfs are a viable model for the companion in the most compact of the X-ray millisecond 
pulsar systems. 

In this Letter, we report the discovery of SWIFT J1756.9— 2508 (hereafter referred to 
as J1756). with a frequency of 182.07 Hz. It was discovered by the Swift Gamma-Ray Burst 
Explorer Burst Alert Telescope (BAT) as part of the hard X-ray transient monitor program. 
In § 2 we discuss the observations made with Swift and the Rossi X-ray Timing Explorer 
(RXTE). In § 3 and § 4 we discuss the timing and spectral analysis and in § 5 we derive 
limits on the companion mass and compare this source to the other known accreting X-ray 
millisecond pulsars. All errors are at the 90%-confidence level unless otherwise stated. 
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2. Observations and data reduction 
2.1. Swift 

SWIFT J1T56.9— 2508 was discovered using the Swift BAT hard X-ray transient monitor 
(Krimm et al. 2006). The BAT is a large field of view instrument that continually monitors 
the sky to provide the gamma-ray burst trigger for Swift. On average more than 70% of 
the sky is observed on a daily basis. Sky images are processed in near real-time to detect 
astrophysical sources in the 15-50 keV energy band 12 . 

J1756 was detected in the daily full-sky mosaics that are produced in the transient 
monitor. It first rose to detectability (Fig. 1) on 2007 June 7 (MJD 54258) with an average 
rate in Swift BAT of 0.0083 ±0.0015 ct s" 1 cm -2 , corresponding to approximately 37 mCrab. 
The discovery was announced on June 13 (Krimm et al. 2007). The BAT data for all pointings 
when the source was fully coded (visible to all BAT detectors) were processed using the Swift 
analysis tools, which included cleaning out the diffuse background and flux from other bright 
sources in the field of view. 

The discovery of J1756 triggered a target of opportunity (ToO) observation in the Swift 
X-ray Telescope (XRT) beginning at 15:25:32 UT, 2007 June 13. It was intermittently ob- 
served over 13 days for a total exposure of 60 ks. The XRT data were processed with standard 
procedures (xrtpipeline v0.10.6 within FTOOLS in the Heasoft package v.6.1.2), filtering, 
screening, and grade selection criteria (Burrows et al. 2005). The photon counting (PC) data 
collected during the first three segments were corrected for pile-up and for the presence of 
single-reflection rings due to the nearby bright source GX 5—1. 13 For the windowed timing 
data we extracted source events in a square region with a side of 20 pixels. Ancillary response 
files were generated with xrtmkarf and account for different extraction regions, vignetting 
and Point-Spread Function (PSF) corrections. The light curve (Fig. 1) was corrected for 
PSF losses and for vignetting. The best fit XRT position is shown in Table 1. 


2.2. RXTE 

The discovery of J1756 also triggered a ToO observation in the Rossi X-ray Timing 
Explorer (RXTE), starting at 18:42:00 UT, 2007 June 13. Because GX 5—1 is within the 
~ 1° field of view of the RXTE Proportional Counter Array (PCA) collimator, the PCA 


12 See http:/ / swift.Ksfc.nasa.gov/docs / swift /resnlts/ transients/ for details on the transient monitor. 
1,! XRT observations after MJD ■ r >4'27S air not usable due to severe contamination from GX 5—1. 
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pointing direction was offset by about 18'. primarily eastward of the true pulsar position. The 
PC A light curve (Fig. 1) has been corrected to effective on-axis values using the collimator 
response function, and spectral response functions account for the offset as well. The offset 
pointing direction is within 3' of the Galactic ridge, and so contamination from diffuse 
emission is present. The amount of contamination was found by examining other nearby 
PCA observations, when J1756 was quiescent, and fitting with a simple phenomenological 
model (and assuming it is constant). The contaminating flux was about. ~3 rnCrab. The 
total exposure was «104 ks in two observing programs. 

Archived monitor data from RXTE and Swift were searched for earlier outbursts of 
J1756. The PCA performs regular monitoring observations of the Galactic bulge (Swank & 
Markwardt 2001) and a search of archived PCA bulge scan data shows that the source has 
not been detected in the 2-10 keV band over the past 8.4 years, with an upper limit of about 
5 rnCrab (2-10 keV). 14 This limit is likely to be systematics driven, based on the nearby 
bright source GX 5—1. Archived daily BAT monitor data (15-50 keV) were searched back 
2.5 years and the RXTE All-Sky Monitor light curves (2-10 keV) were searched back 11.4 
years with no evidence of the source in either instrument with 3<r upper limits of « 15 rnCrab 
(BAT) and 50 rnCrab (ASM). Due to observing constraints (mostly for sun avoidance), 
the overall good coverage of the three detectors is « 87%. With these gaps and the fact 
that even the 2007 outburst was not seen strongly in the ASM, it is possible that an earlier 
outburst was missed. Thus we set a tentative recurrence time scale of > 10 yr. 


2.3. Other Observations 

The discovery announcement led to several multi wavelength observations of J1756. A 
possible fading near-IR counterpart was announced with K s = 19.7 on MJD 54270 and K s 
= 21.0 on MJD 54282 (Burderi et al 2007). However an observation with the Westerbork 
Synthesis Radio Telescope during the X-ray outburst (1.4 GHz; MJD 54266.103; Hessels & 
Stappers 2007) showed no sign of coherent radio pulsations. Furthermore. J1756 was not 
detected in observations after the outburst in either 8.7 GHz radio (MJD 54277.371. Possenti 
et al 2007) or 0. 3-8.0 keV X ray (Chandra upper limit of 1.4 x 10 -13 erg cm -2 s -1 . starting 
MJD 54287.059. Papitto et al 2007) observations. 


14 The 2007 outburst of .11750 was not seen in the PCA bulge scan because RXTE spacecraft slewing 
constraints prevented observations in early Juno. The most recent observation before the outburst was 2007 
Juno (j. 
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3. Timing Analysis 

Fig. 1 shows that the outburst of J1756 had a very rapid rise (< 1 d). After this the 
BAT flux was nearly constant for the next several days. For the first several days of XRT 
and PCA observations (up to MJD « 54268.6), the flux in both instruments showed a slow 
decline, which can be fit to a nearly exponential decay, e~ t ^ T , where r = 7.3 ± 0.4 d. This 
matches fairly closely the decay in the BAT light curve. After this date the flux in all three 
instruments began a steeper decline (r = 0.6 ± 0.3 d) and was undetectable even in the PCA 
after MJD 54271, a mere 13 days after the start of the outburst. This light curve is similar 
in form to those of three other X-ray millisecond pulsars: SAX J1808.4— 3658 (Bildsten & 
Chakrabarty 2001) and IGR J00291+5934 (Gallow r ay et al. 2005), and in particular XTE 
J1751— 305 (Markwardt et al. 2002), which showed very similar values of r (Gierlinski & 
Putanen 2005). The three other objects have had more than one outburst, on recurrence 
time scales of ~ 2 — 3 yr, and all outbursts of XTE J1751— 305 and IGR J00291+5934 and 
most outbursts of SAX J1808.4— 3658 have been short (< 15 days), although notably the 2000 
(Wijnands et al. 2001) and 2002 (Markwardt et al. 2005) outbursts of SAX J1808.4— 3658 
were longer and more erratic. Given the non-detection of earlier outbursts, J1756 is alone 
among the fast-decaying millisecond pulsar systems in not showing recurrence. 

We searched for high frequency variability by constructing an FFT-based pow r er spec- 
trum from the PCA data. The spectrum show r ed a significant, narrow excess near 182 Hz 
(Markwardt et al. 2007), which was also found in following observations. The feature w r as 
several mHz wide, which suggested an X-ray pulsar with orbital Doppler broadening. Exami- 
nation of power spectra on short (f»300 s) intervals confirmed a sinusoidal orbital modulation, 
making SWIFT J1756.9— 2508 the eighth known accreting millisecond X-ray pulsar, and the 
one w r ith the longest rotation period. 

We performed precision timing analysis on the PCA data using the techniques described 
in Markwardt et al. (2002). X-ray event data were recorded in the E_125us_64M_0_ls mode, 
and event times were corrected to the solar system barycenter using the position determined 
from XRT imaging (see Table 1). The Z 2 statistic (or Rayleigh statistic) were computed, and 
the orbital/timing parameters were adjusted to produce the largest Z 2 power of the pulse 
fundamental. This procedure is equivalent to radio or X-ray pulse time-of-arrival fitting 
(Buccheri et al. 1983). but with no binning. After adjustment, we found no long term or 
orbital systematic trends. The binary orbit model was the ELL1 from TEMPO software 
version 11. 010. 15 which is designed for nearly circular orbits. The orbital parameters are 
shown in Table 1, and confirm that this is a very compact and circular orbit. The pulsation 


1 ' r ’]itfp://www.atiif.c.sin).au/rosc'arch/pulsar/tc'iupo/ 
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frequencies from MJD 54264.78 to 54267.07. folded on a trial orbital period clearly establish 
the orbital frequency modulation. The pulsed semi-amplitude is a nearly constant ^6% 
throughout the observations. No pulsed signal was seen in the 15-50 keV band from 250 s of 
BAT event starting starting at MJD 54266.065, with a 2a limit of < 13 %. With the present 
data we can place only weak upper limits on the pulse frequency derivative (/), since the 
time baseline is quite short. 

The data from both the PCA and XRT (which cover 100% of the orbit) were searched 
for modulation at the orbital period. The PCA light curve (2-10 keV) shows no modulation 
at the orbital period, after the long term trend is subtracted. The 2cr upper limit to orbital 
modulations is < 1.2 % from the PCA data. There is also no significant modulation seen in 
the XRT data. Furthermore there are no signs of eclipses, dips or bursts. We searched for 
rapid variability by constructing power spectra of 64 s segments of event data, with a Nyquist 
frequency of 2048 Hz, and averaged contiguous power spectra. While some variability is 
present, especially below ~200 Hz, a detailed timing study is beyond the scope of this 
paper. We do not find any features above 200 Hz which are obvious kiloHertz quasi-periodic 
oscillations (QPOs). For the summed outburst, the 95% upper limits on the fractional r.m.s. 
are about 6% (9%) for 16 Hz (128 Hz) wide QPOs, and a factor of ~2 larger for a single 
observation. 


4. Spectral Analysis 

We fit the Swift/XRT spectrum from each segment in the 0.5-9 keV energy range with 
different models, which included a single absorbed power law (PL; Wilms et al. 2000), and 
an absorbed black body (BB). While a formal f-test cannot distinguish between PL and BB 
models, we favor the PL model because the reduced \ 2 is smaller and closer to unity for the 
PL fit, while the BB fit produces systematic residuals at both low and high energies. The 
derived absorbing column density is always in excess of the Galactic value (N§ — 1.47 x 10 22 
cur 2 , Dickey & Lockman 1990) and there is no indication of a softening of the spectrum 
with time. Representative spectral fits are shown in Table 2. The source was detected in 
the PCA up to energies ~ 30 keV. and the spectrum is consistent with an absorbed power 
law with absorption fixed at the best-fit value of XRT segment 003. The BAT data from 
14 to 195 keV were fit by a simple power law model. The BAT spectrum is softer than the 
spectra found from fits to the PCA alone or XRT alone, suggestive of a cut-off to the power 
law. However joint fits to the BAT. XRT and PCA spectra are not improved with a cut-off 
power law model. 

Given the uncertainty in the spectrum at high energies, we estimate the integrated 
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outburst fluence to be ~ (4.5 ± 0.8) x 10~ 4 erg cm -2 (1-10000 keV). Using the angular 
proximity of the source to the Galactic center to estimate the distance, we can, following 
Bildsten & Chakrabarty (2001), estimate the time-averaged mass accretion rate, M\ = 
(9.3 ± 1.6) x 10 _13 A/3 yr -1 dgraj~ 4 T' 1 Q 1 , where dg is the distance is units of 8 kpc. mi .4 — 
Mx/ (1.4 Mt ) , and we assume a neutron star radius of 10 km. The parameter Tio is the 
recurrence time in units 10 yr, and only T i0 > 1 is allowed by the observations. 


5. Discussion 

In combination, the measured mass function and binary orbital period Pb constrain the 
donor’s mass-radius ( Mq-Rc ) relation, assuming the donor fills its Roche lobe. For SWIFT 
J1756.9— 2508, Pb = 54.7 minutes, a value near the median of the known X-ray millisecond 
pulsars. If we assume a particular value for M x . we can calculate possible values of A/ c from 
the mass function, f x = (A/ c sin i) 3 / (A/ r + A/ c ) 2 , where i is the binary inclination to 
our line of sight and M x and A/ c are the masses of the pulsar and companion, respectively. 
The lack of eclipses allows us to set a limit i < 85° (Paczynski 1971). For this inclination 
and M x = 1.4 A/r. we have A/ C , m , n = O.OO67A/0. Assuming a uniform distribution of 
possible values of cosi we can set a 95% confidence level (C.L.) upper limit on the mass of 
A/ c < 0.022 M© (for M x — l AM-,). The range of possible values of M c up to this limit 
and corresponding companion radius R c are shown in Fig. 2. If we consider a larger neutron 
star of M x — 2.2 A/-, the limits (95% C.L.) are M c < 0.030 A/r and R c < 0.069 P.r ; . This 
radius is less than the minimum radius for low-mass hydrogen rich brown dwarfs of any 
age (Chabrier et al 2000; Bildsten & Chakrabarty 2001, Figure 2), as expected for a binary 
period of < 80 min (Paczynski fc Sienkiwwicz 1981; Nelson et al. 1986). Thus we restrict 
our attention to white dwarf (WD) donor models. 

The dotted and dashed curves in Fig. 2 show evolutionary sequences of ultracompact 
binaries with He donor models originally developed for a study of the AM CVn (double WD) 
class of ultracompact binaries (Deloye et al. 2007). Along these curves, regions of negative 
slope correspond to adiabatic donor evolution. In regions of positive slope the donor is 
cooling and contracting towards the fully-degenerate He WD Mc-Rc relation (see Deloye 
et al., 2007), essentially traced by the lowest dotted line in Fig. 2. The donor's cooling 
is affected by the (very uncertain) efficiency with which irradiating flux generated by the 
accretion — which always dominates the donor's own flux in these ultracompact binaries — is 
thermalized in the donor's upper atmosphere. The impact of irradiation is to extend the 
donor's adiabatic evolution phase and to slow its rate of cooling afterwards. 

Since calculations similar to the Deloye et al. (2007) He-donor study do not yet exist for 
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C/0 donors, we compare the J1756 constraints to a sample of pure C models of Delove A 
Bildsten (2003. .solid curves in Fig. 2); realistic C/O WDs will lie at somewhat smaller Re- 
These models assume fully-convective interior profiles and do not explicitly treat thermal 
transport. However, the physics that produces the transition from adiabatic evolution to 
cooling in He donors does not depend on composition and C/O donors will also experience 
a cooling phase. This means a substantial fraction of the phase space for “hot” C/O donors 
shown in Fig 2 will likely not be evolutionarily accessible (e.g., the nearly vertical regions 
of the log(T c /K) > 6.4 “isotherms : ’), making a C/O donor in J1756 rather unlikely. Thus, 
J1756 likely harbors a He-dominated donor whose thermal cooling has been at least modestly 
slowed. Irradiation provides a natural mechanism that is able to affect the necessary degree 
of slowing, as discussed above. This interpretation places SWIFT J1756.9— 2508 in the 
company of at least two other X-ray millisecond pulsars with very low mass companions: 
XTE J0929— 314, (Galloway et al. 2002) and XTE J1807— 294 (Falanga et al. 2005); although 
J1756 seems to require modest irradiation while these other two sources are compatible with 
the unirradiated cases. 
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Table 1. Timing Parameters of SWIFT J 1756.9—2508 


Parameter 

Value 

Right ascension, a (J2000.0) a 

17 h 56 m 57^35 

Declination, 5 (J2000.0) a 

25° 06' 27" 8 

Barycentric pulse frequency, f ( , (Hz) 

1 82.065804253(72) b 

Pulsar frequency derivative, |/| (Hz s -1 ) 

< 1 x 10 -12c 

Projected semimajor axis, a., sin? (lt-ms) 

5.942(27) 

Binary orbital period, Pi (s) 

3282.104(83) 

Time of ascending node, Tas C d 

54265.28707(6) 

Orbital eccentricity, e c 

< 0.026 c 

Pulsar mass function, f x (10~ 7 A/.») e 

1.56(3) 

Minimum companion mass, M r (10 -3 A/..,) 

i e 6. 7-9.2 f 

Maximum Power, Z% iax 

2061 


a The estimated position uncertainty is 3'/ 5. Position is from the 
XRT and held fixed as a timing parameter. 

b Uncertainties are la in the last quoted digits. 

C 959c upper limit. 

cl Modified Julian days, referred to barycentric dynamical time. 
c Derived parameter. 

f For neutron star masses of 1.4 — 2.2A /». 


Table 2. Sample spectral fits for SWIFT J1756.9— 2508 



MJD 

Model a 

V H 

r [feT] b 

xLi/ dof 

XRT 

54266.06 

PL 

5.44 ±0.41 

2.08 ±0.14 

1.081/200 

XRT 

54266.06 

BB 

2.82 ± 0.23 

[1.30 ±0.05] 

1.185/200 

PC A 

54264.79 

PL 

5.44 

2.00 ±0.14 

0.878/47 

PCA 

54267.46 

PL 

5.44 

1.95 ±0.08 

0.514/48 

BAT 

54258 

PL 

- 

2.3 ±0.7 

0.20/6 

BAT 

54266 

PL 

- 

2.8 ±0.5 

0.38/6 


a P L = power- law: BB=absorbed black body 

b Numbers in column 5 are T for PL fits and kT (in square brackets) 
for the BB fit. 
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MJD 

Fig. 1. — The combined light curve of SWIFT J1756.9— 2508 from Swift BAT (diamonds), 
Swift XRT (filled circles ), RXTE PCA (open circles ), and Chandra (triangle; Papitto et al 
2007). The BAT points show flux in the 15-50 keV band. Chandra the 0. 3-8.0 keV band, 
and the other points the 2-10 keV band. The PCA data is contaminated by emission from 
the Galactic ridge. To correct for this, a constant baseline of 1.24 x 10~ 10 erg cnf 2 s -1 
was empirically determined and subtracted from the PCA data. Data points with < 2a 
significance are shown as upper limits. For comparison with the archival upper limits quoted 
in §2.2, for the PCA (BAT) 1 mCrab ~ 2.4 (5.2) x 10 -11 erg cm -2 s _1 . 
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Fig. 2. — Companion radius (Re) vs. mass (Me) plane, showing the Roche lobe constraints 
for SWIFT J 1756.9—2508 and. for comparison, three other X-ray millisecond pulsars with 
low-mass companions. The source curves trace possible values of the inclination angle, i. The 
solid black dots indicate i = 60° and the open black dots indicate i = 30°. The other curves 
are various white dwarf models. The dotted and dashed lines show evolutionary sequences 
with He donor models calculated as in Deloye et al. (2007) (see text for full explanation). 
The dotted curves show cases where no accretion flux is reprocessed. The dashed curves show 
the case where 109c of accretion flux, as seen by the donor, is thermalized in the donor's 
atmosphere. The numbers next to these curves indicate log tv. where higher values of the 
central degeneracy parameter y c (defined in Deloye et al. 2007) indicate a more degenerate 
and hence compact donor. Within each set of models, evolution proceeds to lower Me and 
less degenerate donors evolve at larger R c . The solid curves are for a pure C composition 
from the models of Deloye Sz Bildsten (2003). The number next to each solid curve indicates 
the log of the temperature in Kelvin used in the model. 




